Real-World Progression-Free Survival as an Endpoint in Lung Cancer: Replicating Atezolizumab and Docetaxel Arms of the OAK Trial Using Real-World Data.

Evaluating cancer treatments in real-world data (RWD) requires informative endpoints. This study replicated the atezolizumab and docetaxel arms of the OAK trial using RWD and compared progression-free survival (PFS) outcomes derived from abstracted physician's notes in RWD (rwPFS) against PFS outcomes derived from the clinical trial PFS (ctPFS). Atezolizumab and docetaxel arms of the phase III OAK randomized controlled trial (RCT; NCT02008227) were replicated in a US nationwide real-world database using selected OAK inclusion/exclusion criteria and propensity score-based adjustment for baseline prognostic variables. Concordance of outcomes was assessed using Kaplan-Meier medians and hazard ratios (HRs). The RWD cohorts comprised 133 patients on atezolizumab and 479 patients on docetaxel. After adjustment, prognostic variables were balanced between RCT arms and corresponding RWD cohorts. The rwPFS and ctPFS outcomes showed better concordance for docetaxel (2.99 vs. 3.52 months; HR: 0.99, 95% confidence interval (CI): 0.85-1.15) than for atezolizumab (3.71 vs. 2.76 months; HR: 0.8, 95% CI: 0.61-1.02). Excluding events labeled "pseudo-progression" from both RWD and RCT improved concordance for atezolizumab (4.24 vs. 4.14 months; HR: 0.95, 95% CI: 0.70-1.25). These findings were robust across sensitivity analyses. Replicating RCTs using RWD and comparing outcomes can help characterize RWD endpoints. Similarity of results between rwPFS and ctPFS at the cohort level may depend on drug category, highlighting the need for further studies to verify and understand when the corresponding outcomes can be compared, including within the same patient.

A meta-analytic framework to adjust for bias in external control studies.

While randomized controlled trials (RCTs) are the gold standard for estimating treatment effects in medical research, there is increasing use of and interest in using real-world data for drug development. One such use case is the construction of external control arms for evaluation of efficacy in single-arm trials, particularly in cases where randomization is either infeasible or unethical. However, it is well known that treated patients in non-randomized studies may not be comparable to control patients-on either measured or unmeasured variables-and that the underlying population differences between the two groups may result in biased treatment effect estimates as well as increased variability in estimation. To address these challenges for analyses of time-to-event outcomes, we developed a meta-analytic framework that uses historical reference studies to adjust a log hazard ratio estimate in a new external control study for its additional bias and variability. The set of historical studies is formed by constructing external control arms for historical RCTs, and a meta-analysis compares the trial controls to the external control arms. Importantly, a prospective external control study can be performed independently of the meta-analysis using standard causal inference techniques for observational data. We illustrate our approach with a simulation study and an empirical example based on reference studies for advanced non-small cell lung cancer. In our empirical analysis, external control patients had lower survival than trial controls (hazard ratio: 0.907), but our methodology is able to correct for this bias. An implementation of our approach is available in the R package ecmeta.

Replication of Overall Survival, Progression-Free Survival, and Overall Response in Chemotherapy Arms of Non-Small Cell Lung Cancer Trials Using Real-World Data.

PURPOSE: The utility of real-world data (RWD) for use as external controls in drug development is informed by studies that replicate trial control arms for different endpoints. The purpose of this study was to replicate control arms from four non-small cell lung cancer (NSCLC) randomized controlled trials (RCT) to analyze overall survival (OS), progression-free survival (PFS), and overall response rate (ORR) using RWD. PATIENTS AND METHODS: This study used RWD from a nationwide de-identified database and a clinico-genomic database to replicate OS, PFS, and ORR endpoints in the chemotherapy control arms of four first-line NSCLC RCTs evaluating atezolizumab [IMpower150-wild-type (WT), IMpower130-WT, IMpower131, and IMpower132]. Additional objectives were to develop a definition of real-world PFS (rwPFS) and to evaluate the real-world response rate (rwRR) endpoint. RESULTS: Baseline demographic and clinical characteristics were balanced after application of propensity score weighting methods. For rwPFS and OS, RWD external controls were generally similar to their RCT control counterparts. Across all four trials, the hazard ratio (HR) point estimates comparing trial controls with external controls were closer to 1.0 for the PFS endpoint than for the OS endpoint. An exploratory assessment of rwRR in RWD revealed a slight but nonsignificant overestimation of RCT ORR, which was unconfounded by baseline characteristics. CONCLUSIONS: RWD can be used to reasonably replicate the OS and PFS of chemotherapy control arms of first-line NSCLC RCTs. Additional studies can provide greater insight into the utility of RWD in drug development.

Building External Control Arms From Patient-Level Electronic Health Record Data to Replicate the Randomized IMblaze370 Control Arm in Metastatic Colorectal Cancer.

PURPOSE: External control (EC) arms derived from electronic health records (EHRs) can provide appropriate comparison groups when randomized control arms are not feasible, but have not been explored for metastatic colorectal cancer (mCRC) trials. We constructed EC arms from two patient-level EHR-derived databases and evaluated them against the control arm from a phase III, randomized controlled mCRC trial. METHODS: IMblaze370 evaluated atezolizumab with or without cobimetinib versus regorafenib in patients with mCRC. EC arms were constructed from the Flatiron Health (FH) EHR-derived de-identified database and the combined FH/Foundation Medicine Clinico-Genomic Database (CGDB). IMblaze370 eligibility criteria were applied to the EC cohorts. Propensity scores and standardized mortality ratio weighting were used to balance baseline characteristics between the IMblaze370 and EC arms; balance was assessed using standardized mean differences. Kaplan-Meier method estimated median overall survival (OS). Cox proportional hazards models estimated hazard ratios with bootstrapped 95% CIs to compare differences in OS between study arms. RESULTS: The FH EC included 184 patients; the CGDB EC included 108 patients. Most characteristics were well-balanced (standardized mean difference < 0.1) between each EC arm and the IMblaze370 population. Median OS was similar between the IMblaze370 control arm (8.5 months [95% CI, 6.41 to 10.71]) and both EC arms: FH (8.5 months [6.93 to 9.92]) and CGDB (8.8 months [7.85 to 9.92]). OS comparisons between the IMblaze370 experimental arm and the FH EC (hazard ratio, 0.85 [0.64 to 1.14]) and CGDB EC (0.86 [0.65 to 1.18]) yielded similar results as the comparison with the IMblaze370 control arm (1.01 [0.75 to 1.37]). CONCLUSION: EC arms constructed from the FH database and the CGDB closely replicated the control arm from IMblaze370. EHR-derived EC arms can provide meaningful comparators in mCRC trials when recruiting a randomized control arm is not feasible.

Will More of the Same Achieve Malaria Elimination? Results from an Integrated Macroeconomic Epidemiological Demographic Model.

Historic levels of funding have reduced the global burden of malaria in recent years. Questions remain, however, as to whether scaling up interventions, in parallel with economic growth, has made malaria elimination more likely today than previously. The consequences of "trying but failing" to eliminate malaria are also uncertain. Reduced malaria exposure decreases the acquisition of semi-immunity during childhood, a necessary phase of the immunological transition that occurs on the pathway to malaria elimination. During this transitional period, the risk of malaria resurgence increases as proportionately more individuals across all age-groups are less able to manage infections by immune response alone. We developed a robust model that integrates the effects of malaria transmission, demography, and macroeconomics in the context of Plasmodium falciparum malaria within a hyperendemic environment. We analyzed the potential for existing interventions, alongside economic development, to achieve malaria elimination. Simulation results indicate that a 2% increase in future economic growth will increase the US$5.1 billion cumulative economic burden of malaria in Ghana to US$7.2 billion, although increasing regional insecticide-treated net coverage rates by 25% will lower malaria reproduction numbers by just 9%, reduce population-wide morbidity by -0.1%, and reduce prevalence from 54% to 46% by 2034. As scaling up current malaria control tools, combined with economic growth, will be insufficient to interrupt malaria transmission in Ghana, high levels of malaria control should be maintained and investment in research and development should be increased to maintain the gains of the past decade and to minimize the risk of resurgence, as transmission drops.

The duration of chemoprophylaxis against malaria after treatment with artesunate-amodiaquine and artemether-lumefantrine and the effects of pfmdr1 86Y and pfcrt 76T: a meta-analysis of individual patient data.

BACKGROUND: The majority of Plasmodium falciparum malaria cases in Africa are treated with the artemisinin combination therapies artemether-lumefantrine (AL) and artesunate-amodiaquine (AS-AQ), with amodiaquine being also widely used as part of seasonal malaria chemoprevention programs combined with sulfadoxine-pyrimethamine. While artemisinin derivatives have a short half-life, lumefantrine and amodiaquine may give rise to differing durations of post-treatment prophylaxis, an important additional benefit to patients in higher transmission areas. METHODS: We analyzed individual patient data from 8 clinical trials of AL versus AS-AQ in 12 sites in Africa (n = 4214 individuals). The time to PCR-confirmed reinfection after treatment was used to estimate the duration of post-treatment protection, accounting for variation in transmission intensity between settings using hidden semi-Markov models. Accelerated failure-time models were used to identify potential effects of covariates on the time to reinfection. The estimated duration of chemoprophylaxis was then used in a mathematical model of malaria transmission to determine the potential public health impact of each drug when used for first-line treatment. RESULTS: We estimated a mean duration of post-treatment protection of 13.0 days (95% CI 10.7-15.7) for AL and 15.2 days (95% CI 12.8-18.4) for AS-AQ overall. However, the duration varied significantly between trial sites, from 8.7-18.6 days for AL and 10.2-18.7 days for AS-AQ. Significant predictors of time to reinfection in multivariable models were transmission intensity, age, drug, and parasite genotype. Where wild type pfmdr1 and pfcrt parasite genotypes predominated (<=20% 86Y and 76T mutants, respectively), AS-AQ provided ~ 2-fold longer protection than AL. Conversely, at a higher prevalence of 86Y and 76T mutant parasites (> 80%), AL provided up to 1.5-fold longer protection than AS-AQ. Our simulations found that these differences in the duration of protection could alter population-level clinical incidence of malaria by up to 14% in under-5-year-old children when the drugs were used as first-line treatments in areas with high, seasonal transmission. CONCLUSION: Choosing a first-line treatment which provides optimal post-treatment prophylaxis given the local prevalence of resistance-associated markers could make a significant contribution to reducing malaria morbidity.

Modelling pathogen load dynamics to elucidate mechanistic determinants of host-Plasmodium falciparum interactions.

During infection, increasing pathogen load stimulates both protective and harmful aspects of the host response. The dynamics of this interaction are hard to quantify in humans, but doing so could improve understanding of the mechanisms of disease and protection. We sought to model the contributions of the parasite multiplication rate and host response to observed parasite load in individual subjects infected with Plasmodium falciparum malaria, using only data obtained at the time of clinical presentation, and then to identify their mechanistic correlates. We predicted higher parasite multiplication rates and lower host responsiveness in cases of severe malaria, with severe anaemia being more insidious than cerebral malaria. We predicted that parasite-growth inhibition was associated with platelet consumption, lower expression of CXCL10 and type 1 interferon-associated genes, but increased cathepsin G and matrix metallopeptidase 9 expression. We found that cathepsin G and matrix metallopeptidase 9 directly inhibit parasite invasion into erythrocytes. The parasite multiplication rate was associated with host iron availability and higher complement factor H levels, lower expression of gametocyte-associated genes but higher expression of translation-associated genes in the parasite. Our findings demonstrate the potential of using explicit modelling of pathogen load dynamics to deepen understanding of host-pathogen interactions and identify mechanistic correlates of protection.

Modelling the benefits of long-acting or transmission-blocking drugs for reducing Plasmodium falciparum transmission by case management or by mass treatment.

BACKGROUND: Anti-malarial drugs are an important tool for malaria control and elimination. Alongside their direct benefit in the treatment of disease, drug use has a community-level effect, clearing the reservoir of infection and reducing onward transmission of the parasite. Different compounds potentially have different impacts on transmission-with some providing periods of prolonged chemoprophylaxis whilst others have greater transmission-blocking potential. The aim was to quantify the relative benefit of such properties for transmission reduction to inform target product profiles in the drug development process and choice of first-line anti-malarial treatment in different endemic settings. METHODS: A mathematical model of Plasmodium falciparum epidemiology was used to estimate the transmission reduction that can be achieved by using drugs of varying chemoprophylactic (protection for 3, 30 or 60 days) or transmission-blocking activity (blocking 79, 92 or 100% of total onward transmission). Simulations were conducted at low, medium or high transmission intensity (slide-prevalence in 2-10 year olds being 1, 10 or 40%, respectively), with drugs administered either via case management or mass drug administration (MDA). RESULTS: Transmission reductions depend strongly on deployment strategy, treatment coverage and endemicity level. Transmission-blocking was most effective at low endemicity, whereas chemoprophylaxis was most useful at high endemicity levels. Increasing the duration of protection as much as possible was beneficial. Increasing transmission-blocking activity from the level of ACT to a 100% transmission-blocking drug (close to the effect estimated for ACT combined with primaquine) produced moderate impact but was not as effective as increasing the duration of protection in medium-to-high transmission settings (slide prevalence 10-40%). Combining both good transmission-blocking activity (e.g. as achieved by ACT or ACT + primaquine) and a long duration of protection (30 days or more, such as provided by piperaquine or mefloquine) within a drug regimen can substantially increase impact compared with drug regimens with only one of these properties in medium to high transmission areas (slide-prevalence in 2-10 year olds ~10 to 40%). These results applied whether the anti-malarials were used for case management or for MDA. DISCUSSION: These results emphasise the importance of increasing access to treatment for routine case management, and the potential value of choosing first-line anti-malarial treatment policies according to local malaria epidemiology to maximise impact on transmission. There is no indication that the optimal drug choice should differ between delivery via case management or MDA.

Asymptomatic Plasmodium falciparum infections may not be shortened by acquired immunity.

BACKGROUND: The duration of untreated Plasmodium falciparum infections is a defining characteristic of the parasite's biology. It is not clear whether naturally acquired immunity (NAI) can shorten infections, despite the potential implications for malaria control and elimination as well as for basic research. METHODS: Data on the presence of P. falciparum msp2 genotypes in six blood samples collected over one year was analysed, together with four samples collected over 1 week, from a cohort in Navrongo (Ghana). Mathematical models assuming either exponential, Weibull, gamma, or log-normal infection durations were estimated separately for six age-groups. The method allowed for varying clonal acquisition and detection rates. RESULTS: The best fitting (Weibull) mean durations were 124 days (children <5 years old), 179 days (5-9 years), and 70-90 days (>10 years). This non-monotonic age pattern is not suggestive of an infection-clearing effect of NAI since immunity increases with exposure, and thus, age. Age-related differences in innate immunity are a more plausible explanation. 21% of blood-stage infections terminated within 1 week, in stark contrast to months of persistence in infections induced in neuro-syphilis patients (malariatherapy data). Age independence in this percentage raises the possibility that this clearance may result from innate mechanisms or genetic incompatibility between hosts and parasites, rather than from NAI. CONCLUSION: In all ages of hosts a substantial proportion of infections are cleared in the first days or weeks of appearance in the blood, while others persist for many months. Although cumulative exposure and NAI increase with age, this does apparently not translate into an increased rate of termination of infections.

Seasonal changes in the antibody responses against Plasmodium falciparum merozoite surface antigens in areas of differing malaria endemicity in Indonesia.

BACKGROUND: The transmission of malaria in Indonesia is highly heterogeneous spatially and seasonally. Anti-malaria antibody responses can help characterize this variation. In the present study antibody responses to Plasmodium falciparum MSP-1 and AMA-1 were measured to assess the transmission intensity in a hypo-endemic area of Purworejo and a meso-endemic area of Lampung during low and high transmission seasons. METHODS: Filter-paper blood spot samples collected from Purworejo and Lampung by cross-sectional survey during high and low transmission season were stored at -20°C. Indirect ELISA assays were carried out using PfMSP1-19 and PfAMA1 antigens. A positivity threshold was determined by samples from local unexposed individuals, and the differences in seroprevalence, antibody level and correlation between antibody level and age in each site were statistically analysed. RESULTS: Prevalence of antibodies to either PfMSP1-19 or PfAMA1 was higher in Lampung than in Purworejo in both the low (51.3 vs 25.0%) and high transmission season (53.9 vs 37.5%). The magnitude of antibody responses was associated with increasing age in both sites and was higher in Lampung. Age-adjusted seroconversion rates showed an approximately ten-fold difference between Lampung and Purowejo. Two different seroconversion rates were estimated for Lampung suggesting behaviour-related differences in exposure. In both settings antibody responses to PfMSP1-19 were significantly lower in the low season compared to the high season. CONCLUSION: Seasonal changes may be detectable by changes in antibody responses. This is particularly apparent in lower transmission settings and with less immunogenic antigens (in this case PfMSP1-19). Examination of antibody levels rather than seroprevalence is likely to be a more sensitive indicator of changes in transmission. These data suggest that sero-epidemiological analysis may have a role in assessing short-term changes in exposure especially in low or seasonal transmission settings.

Comparison of parasite sequestration in uncomplicated and severe childhood Plasmodium falciparum malaria.

OBJECTIVES: To determine whether sequestration of parasitized red blood cells differs between children with uncomplicated and severe Plasmodium falciparum malaria. METHODS: We quantified circulating-, total- and sequestered-parasite biomass, using a mathematical model based on plasma concentration of P. falciparum histidine rich protein 2, in Gambian children with severe (n = 127) and uncomplicated (n = 169) malaria. RESULTS: Circulating- and total-, but not sequestered-, parasite biomass estimates were significantly greater in children with severe malaria than in those with uncomplicated malaria. Sequestered biomass estimates in children with hyperlactataemia or prostration were similar to those in uncomplicated malaria, whereas sequestered biomass was higher in patients with severe anaemia, and showed a trend to higher values in cerebral malaria and fatal cases. Blood lactate concentration correlated with circulating- and total-, but not sequestered parasite biomass. These findings were robust after controlling for age, prior antimalarial treatment and clonality of infection, and over a realistic range of variation in model parameters. CONCLUSION: Extensive sequestration is not a uniform requirement for severe paediatric malaria. The pathophysiology of hyperlactataemia and prostration appears to be unrelated to sequestered parasite biomass. Different mechanisms may underlie different severe malaria syndromes, and different therapeutic strategies may be required to improve survival.

Measurement of Plasmodium falciparum transmission intensity using serological cohort data from Indonesian schoolchildren.

BACKGROUND: As malaria transmission intensity approaches zero, measuring it becomes progressively more difficult and inefficient because parasite-positive individuals are hard to detect. This situation may arise shortly before achieving local elimination, or during surveillance post-elimination to prevent reintroduction. Antibody responses against the parasite last longer than the infections themselves. This "footprint" of infection may thus be used for assessing transmission intensity. A statistical approach is presented for measuring the seroconversion rate (SCR), a correlate of the force of infection, from individual-level longitudinal data on antibody titres in an area of low Plasmodium falciparum transmission. METHODS: Blood samples were collected from 160 Indonesian schoolchildren every month for six months. Titres of antibodies against AMA-1 and MSP-1(19) antigens of P. falciparum were measured using ELISA. The distribution of antibody titres among seronegative and -positive individuals, respectively, was estimated by comparing the titres from the study data (a mixture of both seropositive and -negative individuals) with titres from a (unexposed) negative control group of Indonesian individuals. Two Markov-Chain models for the transition of individuals between serological states were fitted to individual anti-PfAMA-1 or anti-PfMSP-1 titre time series using Bayesian Markov-Chain-Monte-Carlo (MCMC). This yielded estimates of SCR as well as of the duration of seropositivity. RESULTS: A posterior median SCR of 0.02 (Pf AMA-1) and 0.09 (PfMSP-1) person(-1) year(-1) was estimated, with credible intervals ranging from 1E-4 to 0.2 person(-1) year(-1). This level of transmission intensity is at the lower range of what can reliably be measured with the present study size. A Bayesian test for seroconversion of an individual between two observations is presented and used to identify the subjects who have most likely experienced an infection. Furthermore, the theoretical limits of measuring transmission intensity, and how these depend on duration and size of a study as well as on transmission intensity itself, is illustrated. CONCLUSIONS: This analysis shows that it is possible to measure SCR's from individual-level longitudinal data on antibody titres. In addition, individual seroconversion events can be identified, which can be useful in assessing interruption of transmission. Analyses of further serological datasets using the present method are required to improve and validate it. This includes measurement of the duration of antibody responses, how it depends on host age or cumulative exposure, or on the particular antigen used.

The dynamics of natural Plasmodium falciparum infections.

BACKGROUND: Natural immunity to Plasmodium falciparum has been widely studied, but its effects on parasite dynamics are poorly understood. Acquisition and clearance rates of untreated infections are key elements of the dynamics of malaria, but estimating these parameters is challenging because of frequent super-infection and imperfect detectability of parasites. Consequently, information on effects of host immune status or age on infection dynamics is fragmentary. METHODS: An age-stratified cohort of 347 individuals from Northern Ghana was sampled six times at 2 month intervals. High-throughput capillary electrophoresis was used to genotype the msp-2 locus of all P. falciparum infections detected by PCR. Force of infection (FOI) and duration were estimated for each age group using an immigration-death model that allows for imperfect detection of circulating parasites. RESULTS: Allowing for imperfect detection substantially increased estimates of FOI and duration. Effects of naturally acquired immunity on the FOI and duration would be reflected in age dependence in these indices, but in our cohort data FOI tended to increase with age in children. Persistence of individual parasite clones was characteristic of all age-groups. Duration peaked in 5-9 year old children (average duration 319 days, 95% confidence interval 318;320). CONCLUSIONS: The main age-dependence is on parasite densities, with only small age-variations in the FOI and persistence of infections. This supports the hypothesis that acquired immunity controls transmission mainly by limiting blood-stage parasite densities rather than changing rates of acquisition or clearance of infections.

Estimating Plasmodium falciparum transmission rates in low-endemic settings using a combination of community prevalence and health facility data.

As some malaria control programs shift focus from disease control to transmission reduction, there is a need for transmission data to monitor progress. At lower levels of transmission, it becomes increasingly more difficult to measure precisely, for example through entomological studies. Many programs conduct regular cross sectional parasite prevalence surveys, and have access to malaria treatment data routinely collected by ministries of health, often in health management information systems. However, by themselves, these data are poor measures of transmission. In this paper, we propose an approach for combining annual parasite incidence and treatment data with cross-sectional parasite prevalence and treatment seeking survey data to estimate the incidence of new infections in the human population, also known as the force of infection. The approach is based on extension of a reversible catalytic model. The accuracy of the estimates from this model appears to be highly dependent on levels of detectability and treatment in the community, indicating the importance of information on private sector treatment seeking and access to effective and appropriate treatment.

Force of infection is key to understanding the epidemiology of Plasmodium falciparum malaria in Papua New Guinean children.

Genotyping Plasmodium falciparum parasites in longitudinal studies provides a robust approach to estimating force of infection (FOI) in the presence of superinfections. The molecular parameter (mol)FOI, defined as the number of new P. falciparum clones acquired over time, describes basic malaria epidemiology and is suitable for measuring outcomes of interventions. This study was designed to test whether (mol)FOI influenced the risk of clinical malaria episodes and how far (mol)FOI reflected environmental determinants of transmission, such as seasonality and small-scale geographical variation or effects of insecticide-treated nets (ITNs). Two hundred sixty-four children 1-3 y of age from Papua New Guinea were followed over 16 mo. Individual parasite clones were tracked longitudinally by genotyping. On average, children acquired 5.9 (SD 9.6) new P. falciparum infections per child per y. (mol)FOI showed a pronounced seasonality, was strongly reduced in children using ITNs (incidence rate ratio, 0.49; 95% confidence interval, [0.38, 0.61]), increased with age, and significantly varied within villages (P = 0.001). The acquisition of new parasite clones was the major factor determining the risk of clinical illness (incidence rate ratio, 2.12; 95% confidence interval, [1.93, 2.31]). Adjusting for individual differences in (mol)FOI completely explained spatial variation, age trends, and the effect of ITN use. This study highlights the suitability of (mol)FOI as a measure of individual exposure and its central role in malaria epidemiology. It has substantial advantages over entomological measures in studies of transmission patterns, and could be used in analyses of host variation in susceptibility, in field efficacy trials of novel interventions or vaccines, and for evaluating intervention effects.

The distribution of Plasmodium falciparum infection durations.

OBJECTIVES: The duration of untreated Plasmodium falciparum infections in naturally exposed human populations is of interest for rational planning of malaria control interventions as it is related to the duration of infectivity. The extent of variability in duration is relevant where transmission is seasonal, and for the planning of elimination efforts. Methods for measuring these quantities from genotyping data have been restricted to exponential models of infection survival, as implied by constant clearance rates. Such models have greatly improved the understanding of infection dynamics on a population level but likely misrepresent the within-host dynamics of many pathogens. Conversely, the statistical properties of the distribution of infection durations, and how these are affected by exposure, should contain information on within-host dynamics. METHODS AND RESULTS: We extended existing methods for the analysis of longitudinal genotyping data on P. falciparum infections. Our method simultaneously estimates force of infection, detectability, and the distribution of infection durations. Infection durations are modeled using parametric survival distributions. The method is validated using simulated data, and applied to data from a cohort study in Navrongo, Northern Ghana. Distribution estimates from exponential, Weibull, lognormal, and gamma models are compared with the distribution of durations in malariatherapy data. CONCLUSIONS: The Weibull model fitted the data best. It estimated a shorter mean duration than the exponential model, which gave the worst fit. The distribution estimates appeared positively skewed when compared with the distribution of durations in malariatherapy data, suggesting that a significant proportion of infections is cleared shortly after inoculation. We conclude that malariatherapy data, the most important source of information on P. falciparum within-host dynamics, may not be representative of the actual processes in natural populations, and should be used with care. Further, conclusions from transmission models assuming exponential infection survival may be biased.

Detectability of Plasmodium falciparum clones.

BACKGROUND: In areas of high transmission people often harbour multiple clones of Plasmodium falciparum, but even PCR-based diagnostic methods can only detect a fraction (the detectability, q) of all clones present in a host. Accurate measurements of detectability are desirable since it affects estimates of multiplicity of infection, prevalence, and frequency of breakthrough infections in clinical drug trials. Detectability can be estimated by typing repeated samples from the same host but it has been unclear what should be the time interval between the samples and how the data should be analysed. METHODS: A longitudinal molecular study was conducted in the Kassena-Nankana district in northern Ghana. From each of the 80 participants, four finger prick samples were collected over a period of 8 days, and tested for presence of different Merozoite Surface Protein (msp) 2 genotypes. Implications for estimating q were derived from these data by comparing the fit of statistical models of serial dependence and over-dispersion. RESULTS: The distribution of the frequencies of detection for msp2 genotypes was close to binomial if the time span between consecutive blood samples was at least 7 days. For shorter intervals the probabilities of detection were positively correlated, i.e. the shorter the interval between two blood collections, the more likely the diagnostic results matched for a particular genotype. Estimates of q were rather insensitive to the statistical model fitted. CONCLUSIONS: A simple algorithm based on analysing blood samples collected 7 days apart is justified for generating robust estimates of detectability. The finding of positive correlation of detection probabilities for short time intervals argues against imperfect detection being directly linked to the 48-hour periodicity of P. falciparum. The results suggest that the detectability of a given parasite clone changes over time, at an unknown rate, but fast enough to regard blood samples taken one week apart as statistically independent.

Recombination in HIV and the evolution of drug resistance: for better or for worse?

The rapid evolution of drug resistance remains a major obstacle for HIV therapy. The capacity of the virus for recombination is widely believed to facilitate the evolution of drug resistance. Here, we challenge this intuitive view. We develop a population genetic model of HIV replication that incorporates the processes of mutation, cellular superinfection, and recombination. We show that cellular superinfection increases the abundance of low fitness viruses at the expense of the fittest strains due to the mixing of viral proteins during virion assembly. Moreover, we argue that whether recombination facilitates the evolution of drug resistance depends critically on how resistance mutations interact to determine viral fitness. Contrary to the commonly held belief, we find that, under the most plausible biological assumptions, recombination is expected to slow down the rate of evolution of multi-drug-resistant virus during therapy.